Formation of Atomic Carbon Chains from Graphene Nanoribbons 
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The formation of one-dimensional carbon chains from graphene nanoribbons is investigated using 
ab initio molecular dynamics. We show under what conditions it is possible to obtain a linear atomic 
chain via pulling of the graphene nanoribbons. The presence of dimers composed of two-coordinated 
carbon atoms at the edge of the ribbons is necessary for the formation of the linear chains, otherwise 
there is simply the full rupture of the structure. The presence of Stone- Wales defects close to these 
dimers may lead to the formation of longer chains. The local atomic configuration of the suspended 
atoms indicates the formation of single and triple bonds, which is a characteristic of polyynes. 
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Nanoelectronics continuously searches for low dimen- 
sional systems which can be used either as nanocontacts 
or nanoconductors. Metallic nanowires, for example, are 
widely studied because they can present quantum con- 
ductance and the capacity to produce atomic chains [Tl|2]. 

Carbon based-systems, such as carbon nanotubes[2 
and more recently graphene and its derivatives [4], are 
another class of materials which have attracted strong 
interest. They present interesting mechanical and elec- 
tronic properties with great potential for applications in 
nanodevices. Useful properties, such as stability, flexibil- 
ity, charge carriers linear dispersion and high mobility, 
spin injection with long relaxation times and correlation 
lengths could lead to applications in spintronics. 

The possibility to join in the same system the fea- 
tures of one-dimensional wires and the properties of 
carbon based materials is a very exciting one. The 
electronic and transport properties of one dimensional 
carbon systems have already been studied by some 
groups [51 [71 El [101 [IB. However, the lack of re- 
liable and effective ways to produce ID carbon chains 
have limited the studies of these systems. Recently, two 
groups fT2l [13] employed a technique similar to the one 
used for the fabrication of metallic quantum wires [M] 
to obtain experimentally stable and rigid carbon atomic 
chains, which brought new attention to this subject and 
opened up a new avenue in the investigation of one- 
dimensional carbon-based systems. 

Even though this recent experimental work obtained 
these chains via removal of carbons atoms using the 
electron beam, one can envisage a situation where 
these chains could be obtained by stretching graphene 
nanoribbons [[T5l [T6] . In the present work we theoretically 
address this question. We perform room temperature ab 
initio molecular dynamics (AIMD)[17 to investigate the 
formation of linear atomic carbon chains from graphene 
nanoribbons. We elucidate the mechanism of formation 
of these chains and show under what conditions it is pos- 
sible to form these wires. 




FIG. 1: Five representative geometries along the molecular 
dynamics simulations for structures (a) Nl; (b) N2; and (c) 
N3. In (b) and (c), it is possible to see the transition from a 
graphene ribbon to a single carbon chain. 



We study graphene nanoribbons with a neck, as shown 
in Fig. 1. The electronic structure and forces were 
obtained via ab initio total energy density functional 
theory[18j calculations [19 . We have investigated three 
types of necks (Figjl]). In one of them the neck has only 
three-coordinated carbon atoms (Figjlja)), and this will 
be labeled from now on Nl. In the other two systems 
there are additional carbon dimers at the neck, which 
have two-coordinated carbon atoms. In one of them only 
one dimer is added to one of the sides (Figjljb)), and 
in the other two dimers are symmetrically added to the 
two sides of the neck (Figjljc)). These systems will be 
labeled N2 and N3, respectively. In Fig{l]we present five 
representative geometries along the molecular dynamics 
for each one of these systems (see also movies S1-S3[2Q]). 
As can be seen, Nl did not form any one-dimensional car- 
bon chain whereas N2 and N3 formed one and two carbon 
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chains, respectively. Thus, the important conclusion is 
that it seems necessary to have dimers of two-coordinated 
carbon atoms in order to form one- dimensional chains. 

In more detail, the system Nl evolved in such a way 
that the external carbon bonds in the neck broke al- 
most simultaneously, followed by a very rapid rupture 
of the central remaining bond. For N2, the rupture 
started at the opposite side from the two-coordinated car- 
bons dimer, and as a zipper mechanism the other three- 
coordinated carbon bonds subsequently broke, creating a 
one-dimensional carbon chain. Upon continuous stretch- 
ing this chain eventually broke. Finally, for N3, the rup- 
ture initiated at the middle of the neck, again at the 
three-coordinated carbon bonds. Two one-dimensional 
chains were formed. Upon further stretching, one of the 
sides presented a reconstruction where a hexagon was 
transformed into a pentagon with the increase of the 
length of the chain attached to this side. As a conse- 
quence, the continuous pulling led to the rupture of the 
smaller chain, which upon retraction towards the bulk 
formed another pentagon. Further stretching caused yet 
another atomic rearrangement, with the first pentagon 
being broken and reconstructed next to the second one, 
resulting in the reduction of the length of the chain, which 
upon further pulling eventually broke. Of course the 
particular details of the dynamical evolution might de- 
pend on the particular details of the initial conditions. 
However, the important conclusion that the formation of 
one-dimensional atomic chains depends on the existence 
of two-coordinated carbon dimers, as discussed above, 
holds true independently of the simulation details [25 . 

In order to understand why this is so, we analyze the 
average bond lengths and bond angles in the central neck 
region, as shown in Fig. 2. These averages were obtained 
during the time of MD simulation for each elongation, 
or strain. For the Nl configuration, the external bonds 
(bonds (1) and (3)) increase faster and they are where 
the rupture will occur. This is connected with the bond 
angle variations, where the external angle (0^3) increases 
significantly. Now, for the N2 configuration, the bonds 
(2), (3) and (4) increase as the system is pulled, with 
bond (4) increasing much faster, which will lead to its 
rupture. The important point is that the bond length 
of the two-coordinated carbon dimer instead of increas- 
ing suffers actually a small decrease. Once again, these 
changes are correlated with the angular variations. There 
is a large decrease in the angle Pi and an increase in the 
angle ai. This indicates that as the system is pulled the 
stress increase can be absorbed via angular variations for 
this two-coordinated carbon dimer, whereas the central 
neck bonds do not have the fiexibility to do so because 
of their local constraints. Thus, all these internal neck 
bonds will break first, forming the one dimensional wire. 
For the N3 system a similar analysis holds true, as can 
be deduced from the values of the bond lengths and bond 
angles shown in Fig. 2. 
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FIG. 2: Relevant bond lengths and angles (schematically 
shown in the central column panels), as a function of the 
strain, for the three systems studied (Nl - top panel; N2 - 
middle panel; N3 - bottom panel). We present the curves all 
the way up to the rupture of the first bond. 
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FIG. 3: (a) Average potential energy as a function of strain 
for the N2 system. Each one of the points corresponds to one 
of the steps shown in the inset, (b) External applied force 
required to cause a certain strain, obtained from Fig. 3(a). 



In Fig. 3(a) we show the average potential energy 
(DFT energy) as a function of strain for the N2 system. 
Each one of the points in this graph were obtained by 
averaging over the MD simulation for a given length Lw 
of the supercell, which corresponds to a certain strain 
(see the inset, where the average over each one of the 
steps corresponds to a point in the graph of energy ver- 
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sus strain). From this graph we can obtain the apphed 
external force required to cause a certain strain, which 
is shown in Fig. 3(b). One can see that the behavior 
is close to harmonic all the way up to the formation of 
the wire, at approximately 12% strain. At this point the 
bonds break and the linear chain is formed. There is a 
relaxation of the system, with a sudden decrease of the 
potential energy. Further pulling of the system from this 
point again causes an elastic energy increase, basically 
absorbed at the one-dimensional carbon chain. When a 
force of approximately 8 nN is reached, there is a rup- 
ture of the chain. This number can be compared to the 
much smaller forces required to break gold atomic chains, 
which are of the order of 1.5 nN. 
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FIG. 4: Average interatomic distances between the suspended 
carbon atoms shown in the atomic model, for structure N2. In 
the bottom right panel a contour plot for the total electronic 
charge density is presented (units of e/A^). 

Regarding the bonding at the one-dimensional chain, 
we show in Fig. 4 the average interatomic distances be- 
tween the suspended carbon atoms, for structure N2. As 
can be seen, there is an alternation between long and 
short distances, indicating the formation of single and 
triple bonds, which is a characteristic of polvvnes[26]. 
The single bonds are the ones closer to the bulk. In Fig. 
4 we also present contour plot for the total electronic 
charge density. The results confirm the assignment of 
single-triple-single bonds at the suspended chain. 

We have also investigated the role of defects on the 
chain formation. Considering that the Stone-Wales 
(SW) is a possible defect to be formed under electron 
irradiation [13 , we have performed a similar simulation as 
described above for the structure shown in Fig. |5| This 
geometry still has the two-coordinated carbon dimer at 
one side of the neck, but has a SW defect at the other 
side. The SW defect can be seen as a rotation of one 
of the bonds of the hexagonal lattice, resulting in a local 
modification of four hexagons into two pentagons and two 
heptagons. In Fig. 5 we present eight geometries along 
the MD simulation (see also movie S4[20 ). The corre- 
sponding average potential energy variation as a func- 




FIG. 5: Eight representatives geometries along the stretching 
process for a nanoribbon with a SW defect. The arrow in 
the last geometry marks the bond that breaks. The top left 
graph shows the variation of the average potential energy as 
a function of strain, where the numbers refer to the presented 
geometries, whereas the bottom left graph shows the external 
applied force required to cause a certain strain. In the bot- 
tom right panel a contour plot for the total electronic charge 
density is presented (units of e/A^). 



tion of strain is shown in the top left panel. As can be 
seen, there are sequences of quasi-harmonic elastic en- 
ergy increases followed by large atomic rearrangements. 
This behavior is similar to what has been observed in 
the evolution of suspended metallic nanowires[Tl [2]. Ini- 
tially there was the rupture of the bonds of the heptagon 
positioned in the middle of the ribbon (geometry II) , fol- 
lowed by the total rupture at the opposite side of the 
two-coordinated carbon dimer (geometry III). The po- 
tential energy decreases from the geometry II to III due 
to the retraction of the system, which releases the accu- 
mulated elastic energy. This is the general mechanism 
related to the potential energy drops. This also resulted 
in a suspended one-dimensional carbon chain. However, 
as the simulation evolved, instead of simply observing 
the breaking of the chain, there was a rupture of one of 
the bonds of the hexagon where the chain was attached. 
This hexagon is the neighbor of the remaining pentagon 
of the SW defect. As a result, a larger suspended chain 
was formed (geometry V). Moreover, as the system was 
pulled the chain moved sideways, indicating that depend- 
ing on the local geometry and bonding, the thermal en- 
ergy might be sufficient to cause the motion of the wire 
(geometries VII- VIII). This migration was driven by the 
smaller spatial gap at the final connecting point, which 
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leads to a smaller stress of the wire. Finally, further 
pulling causes the rupture of the chain, once more at one 
of the bonds attached to the bulk. In this case it was the 
bond connected to the remaining heptagon of the SW de- 
fect. For this longer chain we also observed an alternation 
between long and short bonds, indicating the formation 
of single and triple bonds, which is also supported by the 
contour plot of the total electronic charge density shown 
in the bottom right panel of Fig. 5. 

In conclusion, we have analyzed the processes involved 
in the formation of linear carbon chains obtained by 
pulling graphene nanoribbons. We have shown that a 
crucial aspect necessary for the formation of these chains 
is the presence of dimers composed of two-coordinated 
carbon atoms at the edges of the ribbons at the neck po- 
sition, otherwise there is simply the full rupture of the 
structure. The formation of longer chains may be associ- 
ated with the presence of Stone- Wales defects close to the 
carbon dimers. Finally, in all situations where we there 
was the formation of suspended carbon chains, we have 
observed the presence of single and triple bonds, which 
is a characteristic of polyynes. 
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